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f**"*. ' The success of nonviscous hydrodynamics in describing the collective flow properties of bulk low 

p± observables at RHIC has led to the claim that a novel form of strongly coupled Quark Gluon 
Plasma (sQGP) is created in 200 AGeV Au+Au collisions. This success depends strongly, however, 
on the initial conditions assumed in the calculations. In particular, agreement with data is only 
obtained assuming Glauber nuclear reaction plane (participant) geometry. The KLN model of Color 
, Glass Condenstate (CGC) initial conditions require the existence of nonvanishing viscous effects. 

' We develop an improved model of CGC that is more internally consistent between the central and 

diffuse edge regions. The improved model is shown to describe bulk rapidity distributions for a 
wide assortment of system sizes, geometries and energies. The self consistency forces a specific y/s 
dependence of the multiplicity which leads to surprisingly low particle production predictions for 
Pb+Pb collisions at LHC COM energies. These predictions are similar to those made using simple 
linear extrapolations of multiplicity attempted by the PHOBOS collaboration. 
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PACS numbers: 12.38.Mh; 24.85.+p; 25.75.-q 



I. INTRODUCTION 

> \ 

^^O , One of the most exciting discoveries in high energy nucleus-nucleus collisions at the Relativistic Heavy Ion Collider 
(RHIC) is the large magnitude of the elliptic flow parameter V2(p±). t>2 = (cos(2</>)) characterizes the transverse 
azimuthal distribution dN/dp±d(f>. The impact parameter or participant nucleon dependence of (i^) provides a clear 
signature of the strong collective behavior of the created medium [l|, 0, d, H| . The systematics of v-i (with energy, 
nuclear size, and centrality) measured at RHIC is consistent with a description of the created Quark Gluon Plasma 
(QGP) core (with T(x) > T c w 170 MeV) as a non-viscous "perfect" fluid that evolves under Euler hydrodynamic 
equations of motion from an early thermalization time r ~ 0.6 fm [H, Q until the hadronization hypersurface 
(T(x) = T c ) is reached. Below T c the highly viscous hadronic resonance gas "corona" evolves with conventional 
hadronic transport @, Such a hybrid perfect fluid sQGP core plus nonequilibrium hadron resonance gas corona 
description, however, is particularly sensitive to the assumed initial geometric distribution of partons produced in 
nuclear collisions since V2 oc e± = (x 2 — y 2 ) / (x 2 + y 2 ). Here, the averaging () is over the transverse plane where x and 
y are the transverse spatial coordinates and the initial number density is used as the weight. The success of the hybrid 
approach to describe the v% data is contingent on the validity of the widely assumed Glauber model of participant 
nucleon geometry. 

With Glauber initial conditions, the parton interaction cross sections needed to generate the collective behavior 
seen at RHIC are much larger than can be supported via purely perturbative collisional processes [§]. This has 
led to the recent strongly coupled (s)QGP paradigm, which is supported by the success of gauge-string dual Anti- 
deSitter /Conformal Field Theory (AdS/CFT) [lCj models in describing the low shear viscosity implied by RHIC 
data [n|, M, G3- 

Hydrodynamic equations of motion, of course, must be supplemented by specified initial temperature and flow 
velocity fields. Conventionally, these are assumed to be specified on a thermalization proper time r ~ 0.6 fm 
hypersurface 0, @. The agreement with data mentioned previously is achieved using the optical Glauber model 
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to calculate the initial distribution of produced matter in the plane transverse to the beam direction in conjunction 
with Brodsky-Gunion-Kuhn (BGK) [14,1151 dynamics that distribute initial matter in a peculiar "trapezoidal" manner 
in the longitudinal rapidity direction [g, LL6J . 

There exist, however, alternative and theoretically preferable models for predicting the initial conditions of ultra- 
relativistic heavy ion reactions that take into account the physics of gluon saturation at small Bjorken x. Color Glass 
Condensate (CGC) refers to a set of models that aim to calculate QCD interactions at asymptotically large energies (at 
very small Bjorken x <C 1) where gluonic degrees of freedom are dominant and the large occu pat ion numbers of gluons 
provide a classical field description of the physics; the so-called gluon saturation regime [17|, |l8|, [l9|, 20]. Calculations 
of physical interest (e.g. unintegrated gluon distribution functions, total cross sections) can be attempted perturba- 
tively because the characteristic scale at which non linearities become important, Q s , provides a large momentum 
scale relative to Aq CD . 

An implementation of the Kharzeev-Levin-Nardi (KLN) [21( parameterization of the unintegrated gluon distribu- 
tions from CGC is used in @, [1] to show that CGC/KLN initial conditions, in conjunction with ideal hydrodynamics 
in the QGP stage, severely overpredict the low p± data at RHIC. This is explained by the fact that the CGC approach 
generically predicts a significantly larger initial spatial eccentricity of the produced matter than Glauber/BGK mod- 
els [22, 23]. This larger initial spatial eccentricity is efficiently transferred by ideal hydrodynamic motion into a large 
elliptical anisotropy in momentum space; a large V2- Thus, if the CGC initial conditions are the ones that actually 
reign in HIC at RHIC, dissipative effects would have to be included in QGP evolution because hadronic dissipation 
has already been taken into account [8|. Perturbative QCD cross sections could then possibly generate the amount 
of collectivity required at RHIC. 

There remain, however, serious conceptual problems in the original KLN implementation of the CGC. As pointed 
out in [13, [H, [24], [26j], the definition of Q s in the model is non-universal because it is defined to be proportional 
to the Glauber participant density, which implicitly depends on the partner nucleus. This same problem leads to a 
drop in the nuclear saturation momentum swiftly to zero as one goes to the nuclear edge, as opposed to the correct 
lower limit of a single nucleon saturation momentum [22I l23l | . There are also further problems stemming from the 
uncertainty in the evolution and proper treatment of the saturation momentum as one increases Bjorken x fit). The 
detailed properties of the local Q s as one progresses to the nuclear edge are crucial to the determination of the correct 
initial state predictions from CGC, especially in peripheral collisions. 

In this paper, we extend the factorized KLN (fKLN) approach introduced in [23| (what we call the extended 
local fKLN (elfKLN)) in order to correctly calculate the local density evolution of Q s and consistently calculate the 
multiplicity of charged particles, dN c h/dn, for multiple systems and multiple collision energies. We take special care 
to account for a consistent evaluation of the model for different system types and sizes and show that this reproduces 
the multiplicities in both A-A and d-A collisions at RHIC with no varying of input parameters. The self consistency 
conditions force a particular COM energy dependence on the multiplicities calculated from the model which lead to 
surprisingly low predictions for Pb+Pb collisions at LHC energies. These predictions are similar to ones made from 
simple linear extrapolations with energy of particle production from the PHOBOS collaboration (2?j] • We do not 
consider the high Bjorken x evolution of Q s as it does not effect determinations of bulk multiplicity and is better 
probed with more differential obervables, e.g. detailed directed flow vi at high p± [26| . 



II. THE SATURATION MOMENTUM AND UNINTEGRATED GLUON DISTRIBUTION FUNCTIONS 

The basic object required to calculate physical observables in the fc^-factorized [2|| picture employed by the gluon 
saturation models is the unintegrated Gluon Distribution Function (uGDF), 4>a(x, kj_). In principle, the uGDFs 
possess a Bjorken x dependence determined by nonlinear evolution equations of the CGC theory [l8|, [l9|, [2(| and their 
k± dependence is fixed by a characteristic saturation momentum, Q s (x). In the McLerran-Venugopalan approach [ItJ 
the gluon distribution is suppressed below the saturation scale 4>a ~ l°g(Q|//c^) compared to the perturbative form 
4>a ~ kj 2 . The parameterization of the KLN model as used in 0] is similar to the following Lorentzian form of 4>a,b- 

The QCD coupling, a s , is regulated at lower scales by imposing a maximum value a max = 0.5. The constant k is a 
parameter set to reproduce dN c h/dr] as seen in experiment. The transverse coordinate dependence is implicit in the 
saturation momentum determined numerically for each nucleus. 

In the KLN model, the saturation scale is assumed to be proportional to the participant density in a nuclear 
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FIG. 1: (Color Online) The figure shows the squared saturation momentum as a function of transverse coordinate for both an 
Au nucleus at Bjorken x = 0.01 and a Pb nucleus at Bjorken x — 0.01^22, ~ 10 -4 . The horizontal lines are the equivalent 
nucleon saturation momenta shown to guide the eye. 



collision. 
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Ql,A( x , x ±; h ) = ■^-ot 8 (Q l l A )xG nuc {x,Q: 1 A )p Pait ,A(x±;h) 



(2) 



In Eq. H C F = ?fi=± and x is the collinear momentum fraction given by kinematics. The participant density is 
defined in terms of the usual Glauber thickness functions Ta/b(x±) as in 0, [HI- The integrated Gluon Distribution 
Function (iGDF) is calculated to be (in the perturbative regime), 

xG nuc (x, Q 2 ) = K log( Q2 ± k2 )x-\l - x) n (3) 

In Eq. 02 the x~ x term accounts for the rapid growth of small Bjorken x gluons while the factor of (1 — x) n was 
introduced in KLN to account qualitatively for the rapid depletion of gluons as x — > 1 outside the small x framework 
of the CGC model. The unintegrated gluon distribution shown in Q] is defined such that ^>(x,'k±) = — 

k\ d g X k f^ ~ Q 2 S 9 g X k f' ) ■ This relationship holds in the leading logarithm DGLAP kinematic regime where k\ Q 2 S . 

The KLN model has the freedom to set the normalization of xG(x, Q 2 ) independently of <p(x, k^_) due to the K factor 
in Eq. [3] in the leading log approximation. This variable is set so as to reproduce the experimentally implied value 
of Q s . We will specify the numerical value of all the parameters after discussing the adjusted KLN type saturation 
model that we use in the paper. Our model differs from both the KLN and fKLN implementations of the CGC in this 
case, as we self consistently include the high x depletion factor of (1 — x) n in both 4>a{x) and the determination of 
Q s . Other models only include this factor in the uGDFs, leading to a relative depletion of particle production in the 
high x regions in our model relative to previous ones. This feature will be essential to our description of asymmetric 
nuclear collisions later in this paper. 

As can be seen in Eq. [21 the saturation momentum denned in the KLN model is not universal as it explicitly depends 
on the partner nucleus through its dependence on pp a rt( x J_) [IH HH HE|- The problem at hand is to modify the 
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constitutive equations of the KLN model such that the saturation momentum depends only on the Glauber thickness 
function T A / B (x±) of the nucleus in question and still reproduces the phenomenologically successful participant 
dependence of the multiplicity that follows naturally from the conventional KLN model [Hj]. 

Just such an adjustment is detailed in 23j where they create the fKLN model. They show that the correct 
event averaging of the uGDF involved in a nuclear collision leads to the replacements Q^(/9p a rt) - » Qs(Ta/pa) and 
4> — ► pa(/> (using the notation found in [23|]) where pa(x±) = (1 — (1 — (?nnTa{x-±_)/A) a ) is the probability of finding 
at least one excited nucleon at the position x_l in nucleus A. The inelastic nucleon-nucleon cross section is set to 
vnn{\J~s = 200, 5500) = 42, 70 mb. Explicitly writing out the changes, 

QLiO^xx) = -p^-a s (Q 2 A )xG nuc (x,Q 2 A ) 



(x,k_L;xx) =Pa{*±) — 77T2 — T72 — T"7t2 — ■ ( 4 ) 



a s (Ql A )kl +Q 



Note that we no longer have the problems occurring from the Q s — > limit that were found in [26J because the A—*l, 
or alternatively the T A —* 0, limit leads to a finite saturation momentum equal to that of a single nucleon. 

27T 2 1 

Qs, P ( x ) = -7^o: s {Q 2 s )xG nuc (x, Q ) (5) 

L F a NN 

We use a value of A = 0.288 as in [2l[ and set the value of the constant K = 0.56 in Eq. [3] to reproduce (Q 2 S A (x = 
0.01)) ~ 2 GeV 2 as in 0,[2l[. The averaging is over the transverse coordinate plane using Ta as weight. 

Fig. [U shows the squared saturation momentum as a function of the transverse coordinate at given value of Bjorken 
scaling variable, x. One can easily see that as one nears the nuclear edge, the nuclear saturation momentum tends 
towards the nucleon saturation momentum from above. An added nice feature of the fKLN approach is that the 
nucleon saturation momentum is set uniquely by the low density limit of the nucleus saturation momentum and is 
not a free parameter as in [2lj . 

III. MULTIPLICITY IN SYMMETRIC COLLISIONS 

We first concentrate on nucleon-nucleon collision and can calculate the multiplicity of charged particles in p-p(p) 
collisions using the formula of Gribov-Levin-Ryskin [28[ . 

dN hi 2 n it? \ r p± 

= C^^T J d2k±M(k± + P±)/2; Ql P (^))M(^x P-0/2; Ql P {x 2 )). (6) 

In Eq.[6l the collinear momentum fractions are given by kinematics as x\^ = P± exp(±y)/y / s. Eq.[6]can be integrated 
over pj_ to get the multiplicity of charged particles in p-p collisions, dN c h / PP / dy . 

The value of n ~ 0.02 is set to reproduce data from the UA5 collaboration [29| . One subtlety in the analysis is 
the r ole p layed by the transformation from the rapidity y to the pseudorapidity r\. The conventional Jacobian is used 
as in |2l| but there is the matter of setting the parameter m g /p± found in the Jacobian. We use the experimentally 
determined average transverse momentum in the collision as a measure for pj_ [30j and find that the best fit is found 
when m g ~ 0.35 GeV. We use the experimentally found average p± [30] in the Jacobian for all collisions at ^fs = 200 
GeV and assume that the ratio m g /p± does not change as one goes to higher energies. The calculation for dN c h/dy 
is shown alongside dN c h/drj in the figured to show the effect of the Jacobian. 

Fig. [2] shows the results of evaluating dN c hf pp / dy for different COM energies. One can see that the multiplicity is 
described well at ^/s = 200 GeV except at very high rapidities. This will be a recurring feature of the calculation as 
large rapidites are in the fragmentation region of the collision and are not well modeled by gluon saturation physics. 
The theory has been fit to the UA5 data [29( and all other calculations are predictions as there are no more parameters 
to tune. The predictions for the p-p collisions at LHC energies can also be seen in Fig. [2] 

Once the parameters of the theory have been specified by reproducing p-p data, we move on to consider symmetric 
nuclear collisions; Au-Au collisions at the max RHIC energy of y/s = 200 AGeV and Pb-Pb collisions at the expected 
LHC energy of ^fs = 5500 AGeV. When we move over to a nuclear collision, we must treat particle production locally 
as a density in the plane transverse to the beam. The parameter values, in particular the normalization of the uGDF, 
is set to produce the correct number of particles in a p-p collision. In order to interpret this as a local transverse 
plane density, we note that the particle production in Eq. is implicitly spread out over a transverse area the size 
of the nucleon-nucleon cross section, unn- We can further understand the smearing of particle production over an 
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FIG. 2: (Color Online) The rapidity and pseudorapidity density of produced charged particles in nucleon-nucleon collisions are 
shown as a function of the rapidity /pseudorapidity. The model is tuned to reproduce data from the UA5 collaboration pffll ] and 
then predictions are made for the charged multiplicity at LHC energies of \fs = 5500, 14000 GeV. 




FIG. 3: (Color Online) Left : Charged multiplicity distribution as a function of rapidity /pseudorapidity for different centrality 
bins in yfs = 200 GeV Au-Au collisions at RHIC. Data includes both statistical and systematic errors and is from the BRAHMS 
Collaboration [3lj]. Right: The figure shows charged particle production per participant pair as a function of the number of 
participants in the collision. Data from BRAHMS [H| and PHOBOS @. 



area of <tnn by noting (as shown in Eqs. [2] and [4]) that all nucleons present in an interaction region of size ctnn 
are considered part of a local density that determines the Q s for particle production. We convert Eq. [5] into a local 
transverse space density, 

dN c h, pp _ 1 dN c h iPP . . 

d 2 X±d 2 p_i_dy cr NN d 2 p±dy ' 



6 




FIG. 4: (Color Online) Left : Charged multiplicity distribution as a function of rapidity /pseudorapidity for different centrality 
bins in y/s = 5500 GeV Pb-Pb collisions at LHC. Right: The figure shows charged particle production per participant pair as 
a function of the number of participants in the collision at LHC. 



The generalization of this to a full nuclear A+B collision is immediate. 
dN c h/AB , r . 1 f , 2 f , 2 2 a s (p 2 j_) 



dy 



'-(V~s) = ^ f d 2 ^ f d 2 P± A^|il f X d 2 k x 0A((kx+Px)/2;Q^ tA (xi))MQ*x-Px)/2;C£ /B (x 2 )). 

<TNNW s ) J J £>F P± J ' 

(8) 

We regard Eq. [S] as the defining expression for our extended local factorized KLN (elfKLN) model of the CGC. It 
includes the consistent local A — ► 1 limit implemented by fKLN in addition to the self consistent high x depletion 
factor of (1 — x) n included in both the uGDFs and Q s determination. This change makes it possible to set a universal 
normalization of the uGDFs such that the adjustments in going from <fi p to 4>a are twofold; we replace Q s ,p with Q s ,a 
and introduce a factor of pa in the nuclear distribution. 

The left panel in Fig. shows the predictions for dN c h I dy and dN c h/dr) in Au-Au collisions at maximum RHIC 
energies along with data from the BRAHMS Collaboration [31j as a comparison. One can see that we achieve a very 
good description of the data over all centrality classes with the exception of the regions near the projectile and target 
rapidity where fragmentation effects dominate. The centrality cuts are imposed via a geometric Glauber analysis as 
in [2l[. We do not take into account the effects of the event by event fluctuation in geometry caused by considering 
the Monte Carlo Glauber analysis which has been shown to be important in explaining system size dependence of 
observables for smaller systems (especially eccentricity and elliptic flow with fluctuations) [32|, [HI, HH • 

The right panel in Fig. [3] shows more clearly the system size dependence of particle productions at RHIC. It shows 
the charged particle production per participant pair as a function of iVp art for a number of different pseudorapidities. 
We get a good description of the particle production, as seen by the agreement between our calculations and data 
from BRAHMS [U and PHOBOS Note that, unlike previous attempts at reaching the N Pa , Tt /2 -> I limit 
there is no problem in this calculation due to the consistent treatment of the local Q s . The calculation extrapolates 
right to the charged multiplicity in a single p-p collision as one reduces -/Vp art . This agreement is non trivial and is 
a good check of our approach as all parameters have already been fit and self consistency demands that we get the 
correct limit. 

The left panel of Fig. [4] shows the rapidity and pseudorapidity densities calculated from the integral of Eq. [8] in the 
case of a Pb-Pb collision at \/s = 5500 AGeV. The centrality cuts are performed using optical Glauber models as 
in [2lT |. The predictions for LHC multiplicity are surprisingly lower than some previous predictions from saturation 
physics [2l|], even the ones using an fKLN type model (25[. The cause of this difference can be tracked to the form 
for Eq. [8] imposed upon us by the self consistency we require between p-p and A-A collisions. Since we have a 
purely local description of particle production (which (2lT | does not include), we are forced to include the factor of 
1/o~nn to maintain consistency between Eqs. [5] and [5] Therefore, we are sensitive to two processes that increase 
total particle production as one increases COM energy; I) the increase of gluon production per unit area caused by 
small x evolution, and 2) the increase of the underlying interaction cross section o~nn{^/s) used as an "averaging 
area" for coherent interaction (as can be seen by the presence of <jnn in the definition of the local density Ta(xj_)). 
Previous calculations include I) but not 2) and are thus typically larger by a factor of 70/42 ~ 1.66 than the current 
predictions. 
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FIG. 5: (Color Online) A schematic figure showing the changes in the local geometric averaging as one goes from Au+Au 
collisions at RHIC energies to Pb+Pb collisions at LHC energies. Both the increase in size of the nucleus and the typical 
coherence area are shown. 

The reason for this difference is schematically shown in Fig. [5] As one increases yfs, one naturally lowers the typical 
x value being probed in the interaction and thus increases the typical saturation scale. This leads to more particle 
production at higher energies as compared to lower energies. The coherence effect, however, is implicitly spread out 
over an area the size of ctnn(\/~s)- As schematically shown in Fig. [5l this increase in the coherence region limits 
the number of coherent production points one can "fit" into the interactions region, thus leading to a suppression in 
overall particle production. 

Our elfKLN model consistently includes both effects mentioned. We systematically isolate the production per unit 
area by using Eq. [5] and thus can consistently explain both p-p and A-A data at RHIC energies, which previous 
approaches are unable to do without invoking an arbitrary enhancement for their p-p predictions [2lL 12a ]. The 
tendency is similar to the one found at RHIC and we note the consistent p-p system limit as we saw in the RHIC 
calculations. The right panel shows the particle production per participant pair at LHC energies as a function of 
A/part ■ 

IV. MULTIPLICITY IN ASYMMETRIC COLLISIONS 

Eq. [8] can be directly generalized to the case of asymmetric collisions by inserting the Glauber thickness functions 
of the different collision partners used. For the particular case of deuteron-Au collisions at maximum RHIC energies, 
one needs to use a model of the nuclear density associated with the deuteron nucleus. We use the Hulthen distribution 
to describe the deuteron, using both s— and d— wave contributions as in [2l[. Once the distribution is determined, 
we calculate the Glauber thickness function Ta(xj_) and proceed as before. 

Fig. shows the comparison of our calculation with data from the PHOBOS Collaboration [38| for a number of 
different centrality cuts as well as for the minimum bias data. The first thing to notice is that there is very little 
agreement in the total number of particles produced at any centrality. The theoretical calculation consistently under 
predicts the data, so much as predicting a multiplicity for a centrality of 80 — 100% that is below the multiplicity for 
a p-p collision. An interesting point to note, however, is that the shape of the overall data is very well reproduced 
just from the asymmetry of saturation momenta in the collision. The presence of the factor of (1 — x) n in the 
determination of Q s is essential to the reproduction of this asymmetry as it causes a precipitous drop in Q s at 
high rapidities. Previous models plj achieve asymmetric results by supplementing their far too slowly varying Q s 
determinations with a leading logarithm evaluation for the p± integral in Eq. [5] This picks the points in phase space 
that possess the most asymmetric values of Qs between the two collision partners. 
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FIG. 6: (Color Online) The Figure shows a comparison of the multiplicity as a function of the rapidity /pseudorapidity for 
different centralities as measured by the PHOBOS Collaboration [38J and calculated in this paper. Dashed lines signify dN c h/dy 
while solid lines are dNch/dr/. 



Centrality 


(^Part) op 


(^art)op 




(^Part)e, 


- 20% 


11.72 


1.99 


13.5 ± 1.0 


2.0 ±0.1 


20 - 40% 


8.21 


1.91 


8.9 ±0.7 


1.9 ±0.1 


40 - 60% 


4.59 


1.63 


5.4 ± 0.6 


1.7±0.2 


60 - 80% 


1.96 


1.00 


2.9 ±0.5 


1.4 ±0.2 


80 - 100% 


0.47 


0.29 


1.6 ±0.4 


1.1 ±0.2 


minbias 


5.37 


1.37 


6.6 ±0.5 


1.7 ± 0.1 



We note that the determination of centr ality cuts in our calculations was performed using an optical Glauber analysis 
as in [21]. As previously shown in [2lT |34| the deviations between a mean field treatment as in the current paper 
and a fully fluctuating Monte Carlo treatment as performed by the experiment is significant, specially when one is 
studying smaller systems. One way to determine the difference between the two approaches is to tabulate the average 
number of participants calculated in each, as shown in the table above. The experimental data in the table is taken 
from the PHOBOS Collaboration (38[. One can see that there is an appreciable difference in the participant number 
between the two approaches, with the fluctuating approach isolating larger systems (i.e. smaller impact parameters) 
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FIG. 7: (Color Online) The Figure shows a comparison of the multiplicity per participant pair as a function of the rapid- 
ity /pseudorapidity for different centralities as measured by the PHOBOS Collaboration [38| and calculated in this paper. 
Dashed lines signify dN c h/dy while solid lines are dN c h/dr]. 



for each centrality cut. The table also explains the peculiarity found in Fig [6] for the 80 — 100% centrality cut. The 
optical Glauber actually gives total participant number of less than 1 for that centrality cut, thus leading to particle 
production below the p-p level. 

Given the different determinations of centrality in the approaches, a good way to compare the data to the calculation 
is to divide out by the overall number of participant pairs. This controls for the difference in the underlying impact 
parameter differences and poses the problem in terms of particle production per participant pair. The analysis is 
presented in Fig. [7] and it can be seen that the comparison is a lot more favorable now that differences due to 
centrality determination have been removed. There is good agreement between data and theory, specially in the mid 
to forward rapidity region. There is still some discrepancy as one moves backward in rapidity towards the nucleus. 
Intranuclear cascading increases the multiplicity of charged particle in the nucleus fragmentation region. Furthermore, 
we have neglected to take into account the geometric scaling region of the uGDF at low x [3|| which has been shown 
to be important in d-Au collisions at RHIC, at least off mid rapidity 40]. Both these effects are larger for more 
central collisions, a fact that is borne out in Fig. [7] where the discrepancy between data and theory is reduced as one 
moves towards more peripheral collisions. 

Having achieved a good description of d-Au data at RHIC we move on to a prediction for p-Pb collisions at LHC 
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FIG. 8: (Color Online) The figure shows the rapidity distribution of charged particles produced in a p-Lb collision at \fs = 5500 
GeV. Dashed lines denote dN c h/dy while solid lines are dN c h/drj. The open squares are minimum bias predictions. 

energies. Eq. [6] easily generalizes to an asymmetric p-Pb collisions by replacing one of the uGDFs with one that 
uses the saturation momentum of a nucleus. Fig. [5J shows predictions of a p-Pb collision for different centralities at 
y/s = 5500 GeV. The overall charged particle production is smaller by a factor of about 2 from p-Lb predictions 
made using the conventional KLN model (2~H |. One must note, however, that this calculation makes centrality cuts 
using an optical Glauber model calculation and thus will run into similar problems to the d-Au calculations when 
being compared to data. The problems can be seen explicitly in the the 80 — 100% centrality bin in Fig. [5J where the 
particle production is at a level lower than p-p at y/s — 5500 GeV (Fig. [2J. The minimum bias prediction is most 
robust, as was the case with the d-Au predictions. We present this prediction as a good estimate of the expected 
multiplicity, one that will be better compared to data when presented as particle production per participant pair. In 
general, a complete description of asymmetric nuclear collisions may best be attempted using a model that can better 
accommodate fully fluctuating geometries, such as the Monte Carlo KLN (MCKLN) model detailed in [23| . 

V. CONCLUSIONS 

The recent claims supporting the sQGP paradigm at RHIC are heavily dependent on the success of models using 
ideal hydrodynamics to model the evolution of the QGP stage of bulk evolution at RHIC 0, H|. This success is 
inherently linked to the choice of initial conditions assumed in the models, and an initial state derived from gluon 
saturation / CGC models does not work well in conjunction with ideal hydrodynamics to explain RHIC data [8j. This 
is due to the larger initial spatial eccentricity produced in gluon saturation models which is naturally converted into 
a larger elliptic flow parameter, V2 [12] • A better understanding of the origin and robustness of this large eccentricity 
is needed. In order to do this we need to understand the A — > 1 and Ta — * limits of the theory, which is equivalent 
to creating a model that works in describing not just central A-A collisions but also p-p and p-A collisions in one 
consistent model. This is especially important as recent work indicates that models that have more realistic nuclear 
edge characteristics lead to a spatial eccentricity that lies somewhere between the pure Glauber/BGK and KLN model 
eccentricities [23l[23|. 

We use the elfKLN model, our extension of the fKLN model in [23], consistently applied to p-p, p-A and d-A 
collisions along with the usual description of A-A collisions. We find that the model does an excellent job consistently 
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describing particle production in symmetric nuclear collisions from p-p all the way up to A-A collisions. We fix 
our parameters to reproduce particle production in p-p collisions and predict the multiplicity in Au-Au collisions at 
RHIC, finding good agreement with the data. The model does very well in treating the nuclear edges consistently, 
as evidence by Fig. where we show a consistency check between our Au-Au and p-p calculations. We also predict 
the multiplicity for Pb-Pb collisions at LHC energies, which are found to be surprisingly small. This is due to the 
increasing "coherence area" of the interaction, which we consistently build into our model. 

We move on to calculate asymmetric nuclear collisions (d-A, p-A) in an effort to correctly model the edges of 
peripheral A-A collisions. We find that absolute agreement with the data is difficult to attain if one uses optical mean 
field methods to obtain the centrality cuts. However, good agreement is achieved with the d-Au RHIC data if one 
accounts for the variation in the underlying geometry by comparing particle production per participant pair. 

In closing, we propose the clfKLN model as a robust model for bulk CGC physics that has the properties of being 
explicitly factorized and well defined in the diffuse nucleus limit. We do, however, recognize that it it might not be 
applicable in highly asymmetric (and thus smaller system size) nuclear collisions, where a fully fluctuating model such 
as the MCKLN of HH might be more appropriate. Our prediction does not include effects such as running coupling 
and pre asymptotic terms in the evolution of the saturation momentum with energy [36L lo^]. These adjustments are 
left for future work. 
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